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Abstract 

We critically discuss the recent observations of the binary system at the center 
of the bipolar planetary nebula Henize 2-428. We find that the proposed 
explanation of two equal-mass degenerate objects with a total mass larger 
than the Chandrasekhar limiting mass that supposedly will merge in less than 
a Hubble time, possibly leading to a SN la, is controversial. This hypothesis 
relies on the assumption that the variability of the He II 5412 A spectral 
line is due to two absorption components. Instead, we propose that it can 
be accounted for by a broad absorption line from the central system on top 
of which there is a narrow emission line from the nebula. This prompted 
us to study if the binary system can be made of a degenerate star and a 
low-mass main sequence companion, or of two degenerate objects of smaller 
mass. We find that although both scenarios can account for the existence of 
two symmetric broad minima in the light curve, the second one agrees better 
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with observations. We thus argue that the claim that Henize 2-428 provides 
observational evidence supporting the double-degenerate scenario for SN la 
is premature. 

Keywords: Planetary nebulae, Stars: AGB and post-AGB, Supernovae, 
White dwarfs 


1. Introduction 


Thermonuclear, or Type la supernovae (SNe la), are the result of the 
explosion of carbon-oxygen white dwarfs. Despite their well known observed 
properties, the nature of the progenitor systems that produce a SNe la event 
has not been hitherto elucidated, and several scenarios have been proposed, 
none of which gives a satisfactory answer to all the abundant observational 



of some of these scenarios. 

As there is no consensus on which are the SN la progenitor(s), it is cru¬ 
cial to refer to all scenarios (or categories of scenarios) when confronting 
them with observations. We list them in alphabetical order, and cite only 
a few references for each scenario: a) Th e core-degenerate (CD) scenario 


(Livio & Riess, 2003: 

Kashi & Soker. 

2011: Soker et ah. 

20131. bl The double- 

degenerate (DD) scenario (e.g., 

Webbink 

1984: 

ben & Tutukov 1984 

1), cl 

The double-detonation (DDetl mechanism (e.g., 

Wooslev & Weaver 

1994; 


nario 

elT 


Livne fe Arnettlll995l: IShen et al. 120131. dl Th e single-degenerate (SDf see 


e.g., Whelan fe Iben 1973 : Nomoto 19821: Han fe Podsiadlowski 2004 1 


2009 


re recently proposed singly-evolved star (SES) scenario ( Chio si et al. 


2015 ) , and f) The WD-WD collision (WWQ scenario (e.g., 


Thompson! 20 111: Kushnir et ah 2013 : Aznar-Siguan et al. 


Raskin ct al. 


201311 . 


Since all these scenarios involve white dwarfs, all progenitors evolve through 
one or two planetary nebula (PN) phases. Accordingly, one of the pieces of 
evidence that would help in constraining SN la scenarios is to study PNe. 
Furthermore, in some cases SN la have been even claimed to take place in¬ 


side p lanetary nebulae (e.g., iDickel fe Joneslll985l: iTsebrenko fe Sokerl 12013 


2015|), a process ter med SNIP. 


In a recent paper ISantander-Garcfa et al.l (1201511 analyzed the central bi 


nary system of the planetary nebula Henize 2-428 (iRodrfguez et al.l. 2R01 
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Santander-Garcia et all 1201 if ). ISantander-Garcfa et al.l (120151 ) found that 
the light curve of this PN shows two nearly identical broad minima, indicat¬ 
ing significant tidal distortion of the components of binary system, and that 
there is an absorption line of He II 5412 A that varies with time. Given 
that the two minima of the light curve are practically identical, they assumed 
that they are caused by a binary system composed of two equal-mass stars 
of the same type, and found the temperature, radius, and luminosity, of the 
two stars to be almost identical. They further argued that most likely these 
are two degenerate stars, i.e., white dwarfs or cores of post-asymptotic giant 
branch (AGB) stars, on their way to become CO white dwarf s . As t he com¬ 
bined mass in this model is 1.76 M 0 , ISantander-Garcfa et al.l (120151 ) further 
argued that these two stars will merge to form a SN la in the frame of the 
DD scenario. 


Here we critically discuss the explanation of lSantander-Garcfa et al.l (120151) 
As we ex plain in Sect. [2] we find the interpretation of the observations of 
Santander-Garcia, et al.l (120151 ) to be plausible, albeit other possibilities are 
conceivable. In Sect. H we relax the assumptions of these authors and we 
propose alternative models of the binary system. The first of these models 
consists of binary system in which only one of the components is a degener¬ 
ate star, while the secondary star is normal non-evolved star. The second of 
the models involves two non-identical degenerate stars, but with a combined 
mass smaller than the Chandrasekhar limiting mass. A short summary is 
given in Sect. [1 


2. Preliminary considerations 

2.1. A binary system made of two identical stars? 


Santander-Garcia et al.l (120151) argue for a binary system composed of 
two stars having the same mass, 0.88 ± 0.13 M©, the same luminosity, ~ 
420 Lq at a distance of 1.4 kpc, and the same radius, 0.68 ± 0.04 R Q . This 
implies that the two stars are at the same evolutionary stage. However, 
any small difference in the main sequence mass will turn to a large one on 
the asymptotic giant branch (AGB). An AGB star having a core of O.88M 0 
burns hydrogen at a rate of ~ 2 x 10“' M 0 yr” 1 (e.g., Paczvnskil 1970 ). 
For a difference in mass between the two cores < 0.02 M 0 the difference of 
evolutionary times between the post-AGB stars should be ^ 10 5 yr. This 
requires a mass difference on the main sequence of AM /M ^ 10“ 3 , depending 
on the initial mass of the stars. 
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It could be argued that there are other bina ry systems w ith almost iden¬ 
tical c omponents, known as twin bi naries f Lucv fe Ricco . 19791) . Specif¬ 
ically, Pinson ne ault & Stanekl (12006!) studied 21 detached eclipsing bina¬ 
ries in the Small Magellanic Cloud and found that 50% of detac hed bina¬ 
ries have companions wit h very similar masses. However, Lucvl ( 2006 1 and 
Cantrell fe Douganl ( 2014!) concluded that there is a strong observational bias 
that affects spectroscopically selected binary stars, and that the apparent 
overabundance of twin binaries does not reflect their true population. In 
summary, the case for a twin binary is possible, but unlikely, hence moti¬ 
vates us for a careful reexamination of such a claim. 


2.2. Stellar properties 

As mentioned, in the model proposed by ISantander-Garcia et ahl (1201511 
each star is a post-AGB star with a mass of 0.88 M 0 . When a post-AGB of 
that mass fades to a luminosity of ~ 1O 3 L 0 its radius is already ~ 0.02 R® 
(e.g., iBloecker fe Schoenbernerlll99lll. This radius is about 3 0 times smaller 
than the radius suggested by iSantander-Garcia et al. (120151 ). This poses a 


serious problem to their model. 

In the first of our models we investigate a case where the luminosity of 
the system is due to just one star, and the luminosity of the companion is 
negligible — see below for more de t ails. At a distance of D — 1.4 kpc as 
deduced bv ISantander-Garcia et al.l (1201511 the luminosity is ~ 850 L Q . This 
can be a star whose evolution was truncated on the upper red giant branch 
(RGB), when its core mass was only M\ & 0.45M 0 , or on the lower AGB 
when its core mass was ~ 0.5 M 0 . On the other hand, if the distance is 
larger, say D = 1.8 kpc, the luminosity is ~ 1.4 x 1O 3 L 0 . This can be a 
star whose evolution was truncated on the lower AGB, when its core mass 
was only Mi 0.52 — 0.55 M 0 . In our proposed model the companion that 
terminated the RGB or the AGB evolution of the primary component is a 
main sequence star of ~ 0.3 — 0.5 M 0 . In the second of our models we assume 
that indeed both stars are post-AGB stars but we allow the stars to have 
different physical parameters, namely different masses, effective temperatures 
and luminosities. 
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while the most recent determination of lFrew et al.1 (120151 1. using the Ha sur¬ 
face brightness-radius relation is also 2.7 kpc. Based on these values we will 
scale our expressions with two distances, D = 1.4 kpc and D — 1.8 kpc, as 
the value adopted by Santander-Garcfa et ah ( 2Q15h was obtained from their 
fit to the properties of the binary system, which is questioned here. 

2.3. Light curves and spectrum 


The arguments of ISantander-Garcia et al.1 (120151 1 for their claim of a bi¬ 
nary system of equal-mass stars at the same evolutionary stage are the nearly 
identical minima in the light curve, and the line profile of the He II 5412 A spec¬ 
tral feature — see their figures 2 and 3. The nearly identical minima of the 
light curve have been suggested to be indicative that both mem bers of the bi 


nary s ystem have very similar masses. Additionally, ISantander-Garcia et al 


(120151 ) found that the He II 5412 A spectral line of Henize 2-428 is vari¬ 
able. They attributed the variability of this line to Doppler shifts of two 
equal-mass stars, and then used two Gaussian absorption profiles to model 
the variation. Consequently, in their joint analysis of the light curve and the 
spectrum they fo rced the mass ratio q = of the binary system to be 

1. Furthermore, Santander-Garcfa et al.l ( 2015 1 did not model the spectra 
of both components of the binary system, since they were not able to mea¬ 
sure surface gravities for each one of the individual binary members. Finally, 
they assumed that both stars are at the same evolutionary stage. All of these 
assumptions are critical in their analysis. 

In particular, it must be stressed that even if the mass ratio is close to 1, 
the nature of the stars can be very different, and that the lack of determina¬ 
tions of surface gravities leaves room for alternative explanations. In particu¬ 
lar, a close look at figures 2 and 3 of ISantander-Garcia, et al.l (1201511 suggests 
that the spectrum of Henize 2-428 can be explained by assuming that there is 
an emission line on top of the absorption profile. We therefore examine a pos¬ 
sible alternative interpretation where the line profile is the result of a wide 
absorption line belonging to the primary star, and a narrow em ission line 
coming from the compact dense nebula reported by lRodrfguez et al. (120011 1. 
or which originates even much closer to the star from the wind itself. In 
this alternative explanation both emission and absorption lines change with 
orbital phase. This is not unusual. Many cent ral stars of planetary nebulae 


show He II absorption lines (e.g., Weidmann fe Gamen 201lh . The emission 


l ine is seen in som e nebulae, e.g., in Abell 48 whose ce ntral star is a WN star 
( Todt et all . 2013 1. Most interestingly, in the study of Weidmann fe Gamen 
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(201l|) there are several PNe that show a wide He II 5412 A absorption 
line with a weak emission feature in the center of the wide absorption line. 
This forms a spectral structure similar to that of Henize 2-428. The most 
noticeable examples of this are the PNe He 2-105, He 2-434, and to some 
degree SP 3 and PC 12. All these central stars are O stars. To these we 
add the PN Pa 5, for which the central star shows a wide He II 5412 A ab¬ 
sorption line and the res olved nebula shows a narrow He II 5412 A emission 
line (IGarcia-Diaz et all 2 0141 ). Another relevant case is the central star of 
the Eskimo planetary nebula (NGC 2392). The He II 5412 Aline s hows a 
structure similar to that of Henize 2-428 flPrinia fe Urbaneial. I2014T) . This 


struct ure varies on time scales down to about an hour. Prinia fe Urbaneial 


(2014) attribute the variation in the line to a variable distribution of clumps 
in the wind, variations in the velocity held, and/or the mass loss rate, rather 
to a binary star. 

It could be argued that if our interpretation is correct, and the He II 5412 A 
sion line is nebular, a He II 4686 A emission line should be also observed, but 
it is not. However, other PNe, like He 2-105, PC 12, He 2-434 and SP 3 have 
spectra with weak He II 5412 A bumps (which we interpret as an emission 
line) in the center of a broad absorptio n line, whereas the He II 4686 A line 
is in absorption (jWei d mann fe Gamen . 1201111 . On the other hand, in the PN 
M 1-14 the He II 5412 A line is in abso rption but the He II 4686 A occurs 
in emission ( Weidmann & Gamenl . 20 ll! ). These cases show that there is an 
interplay between emission and absorption in the He II lines. This suggests 
that the emission region, like the absorption one, is located very close to the 
star, but not at the same place. In all these cases both the absorption and 
emission features vary with orbital phase. In conclusion, based on the cur¬ 
rent observations the existence of a He II emission line on top of the broad 
absorption featur e cannot be discarded, and nee ds to be considered. 

Furthermore, iSantander-Garcia, et al.l 02015 1 established an upper limit 
on the effective temperature of the members of the binary system of 
40,000 K, based on the absence of a He II 5412 A emission line, but 
there are many PNe whose central stars have effective temperatures larger 
than 40,000 K and do not have He II emission lines. Since we assume 
th at the emission line is indeed present in the spectrum, the upper limit 
of ISantander-Garcia et al.l 020151) to the effective temperature is no longer 


valid, and the effective temperature could be sufficiently large to allow emis¬ 
sion. Consequently, in the following we will scale quantities with an effective 
temperature T eS = 45, 000 K. 


cmis- 
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3. Possible binary models 


According to all the considerations put forward in the previous section, 
here we explore two models in which the mass ratio is ~ 1, but with different 
characteristics. The first of these models involves a semi-detached system in 
which the members of the binary system are not at the same evolutionary 
stage, while for the second mo del we adopt an over contact binary system 
made of two post-AGB stars, as Santander-Garcfa et al. ( 20151) did, but with 
less restrictive assumptions. 


3.1. A semi-detached binary system 

We first examine a binary model where the secondary is a main sequence 
star of mass M 2 ~ 0.3 —O.5M 0 . If we take the primary to be a post-RGB star 
of mass Mi ~ O.45M 0 or a post-AGB of mass Mi ~ O.55M 0 (e.g., iBloecker 
19931 ). then these two cases span a mass ratio of q — My,/Mi 0.6 — 1.0. 
Using the expression for the Roche lobe radius r l from Eggletpn (119831 ). we 
End for the primary star r^i/a = 0.42 and 0.38 for q = M 2 /Mi = 0.6 and 
1.0, respectively, where a is the orbital separation. For the secondary star 
we find r L2 /a = 0.34 and 0.38 for q = M 2 /Mi = 0.6 and 1, respectively. For 
an orbital period of P = 4.2 h the orbital separation (for a circular orbit) is 
a = 1.27(M/0.9 M 0 ) 1//3 i? 0 , and the Roche lobe of the primary star is 


T LI = 0.51 


M 


0.9 M 0 


1/3 


r l/ a 
0.4 


Rq, 


( 1 ) 


where M = Mi + M 2 is the total binary mass. From the primary luminosity, 
Li = 845(_D/1.4kpc) 2 L 0 , and for the effective temperature assumed here, 
the primary radius is 


Ri = 0.48 


D 


T\ 


1.4 kpc J \4.5 x 10 4 K 


-2 


Rq. 


( 2 ) 


For a distance of D = 1.8 kpc the primary radius with that effective temper¬ 
ature is Ri = 0.62 R q , but taking a temperature 2) = 50,000 K will make 
the primary just hlling in its Roche lobe. In our model the primary is close 
to hlling its Roche lobe (assuming synchronization). For the secondary the 
Roche lobe radius is in the range rL 2 — 0.42 — 0.5 Rq. Thus, the secondary 
also touches its Roche lobe. 

We consider two possible evolutionary scenarios. We prefer the post-RGB 
one, but cannot rule out the post-AGB one. We first describe the relevant 
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evolutionary tracks on the Hertzsprung-Russell diagram, and then turn to 
discuss the two possible scenarios. The last step of our analysis consists of 
computing the light curve of the proposed binary system. 


3.1.1. Evolutionary tracks 

In the top panel of Fig. [Owe present some relevant evolutionary tracks 
on the Hertzsprung-Russell diagram, together with the approximate location 
of the central star of Henize 2-428 at two distances, D = 1.4 kpc (square) 


and 1 .8 kpc (asterisk), and the position reported bv ISantander-Garcfa et al. 


(12015j ). circle The sequen ces of masses Mi = 0.57 and 0.87 M 0 were taken 


from iRenedo et al.1 (120101) . and are the result of the full evolution of zero- 
age main sequence (ZAMS) stars with masses M z amSi = 1-5 and 5.0 M e , 
respectively, of Solar metallicity. The evolution is followed through all the 
relevant stages, including the hydrogen and helium core burning phases, the 
thermally-pulsing AGB phase, and the post-AGB evolution to the w hite 


dwarf stage. The Mi = 0.504 M 0 sequence was also taken from IRenedo et a l 


(120101 ). and corresponds to a progenitor of mass M Z amsi = 0.85 M 0 with 
Z = 0.001. In addition, the post-RGB evolution of a 0.432 M 0 helium-core 
low-mass white dwarf is included. This sequence is the result of the non¬ 
conservative binary evolution of a star of mass M Z amSi = 1-0 M^ that aban¬ 
dons the RGB before the onset of core helium burning (lAlthaus et all 1201 3 1 ). 
Finally, the evolutionary track of a Mi = 0.84 M 0 post-AGB remnant of a 
A^zamSi = 3.0 M q progenitor with metallicity Z = 0.001 that experiences a 
late thermal pulse (LTP) when its effective temperature was T e g = 10, 000 K 
is shown as well. As a result of the LTP, the post-AGB remnant experiences a 
fast evolution to the blue. After reaching a maximum effective temperature, 
this remnant evolves back to the domain of giant stars. 

Note that the luminosity inferred for Henize 2-428, for the effective tem¬ 
peratures used by ISantander-Garcfa et al.l (120151 ) and by us, is substantially 
smaller than that predicted by post-AGB evolutionary models for masses of 
A/zamSi ^ O.55M 0 . The claimed masses for the two central stars of Henize 2- 
42 8 do not fit the other claimed p roperties of the stars in the model proposed 
by Santander-Garcfa et al.l ( 20151 ). However, the Mi = 0.504 M 0 post-AGB 
star presented here (upper panel) experiences an LTP at T eff = 48, 000 K, 
which brings the remnant back rapidly to the red giant domain and finally to 
the white dwarf stage. Such an evolutionary track covers the general region 
of Henize 2-428 on the Hertzsprung-Russell diagram. 

Finally, the bottom panel of Fig. Q] displays the evolution of a star with 



























5.5 5 4.5 4 3.5 


LogT 


eff 


Figure 1: Relevant evolutionary tracks (see text for details). Marked are the location of 
each of the two identical stars of Henize 2-428 according to the model of Santander-Garcfa 
et al. (2015), filled circle, and the central star in our proposed scenario when the distance 
is D — 1.4 kpc (filled square) and when D = 1.8 kpc (asterisk). 
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A/zamSi = 2.5 M 0 , and Z = 0.01. The evolution of this star was truncated 
on the lower AGB, when its luminosity is L — 1.4 x 1O 3 L 0 , the observed 
luminosity of the binary system for the distance we adopt. At this point mass 
was removed to mimic a common-envelope episode. The final mass of the 
remnant after this intense episode of mass loss is 0.50 Mg. The evolution in 
the Hertzsprung-Russell diagram of this sequence shows several blue loops, 
which are due to successive ignitions of the hydrogen shell. As can be seen, 
this track is also able to reproduce the observed position of Henize 2-428, for 
our adopted distance. 

In conclusion, post-RGB helium-core remnants with stellar masses of 
Mi ~ 0.45 M 0 and low-mass cores of AGB stars that truncate the lower 
AGB with masses of Mi pa 0.50 — 0.55 M 0 can account for the central star 
of Henize 2-428. 


3.1.2. Post-RGB evolution 


Guerrer o et ah (120001 ) studied 15 bipolar PNe, and only in Henize 2- 


428 they detect no H 2 emission, and it was the only PN in their sample in 
which a bright central star was found. It seems there is something strange 
in the evolution of th is PN, w hich here is attributed to the system being a 
post-RGB star. Rodriguez et ah 1120011 ) found the low abundances of most 
elements they study (relative to hydrogen) of Henize 2-428 to be similar to 
those found for PNe belonging to the Galactic halo. They suggested that 
this points to a low-mass progenitor. It is quite possible that the central star 
of Henize 2-428 is a post-RGB star orbited by a low-mass main sequence 
star. The post-RGB scenario is compatible with the low nebular mass, as 
the initial stellar mass in this scenario is MzamSi ~ 1 M 0 , and the post-RGB 
mass is Mi pa 0.45 M 0 . 


3.1.3. Post-AGB evolution 

In this proposed alternative scenario the evolution of the primary star 
is truncated on the lower AGB when its luminosity is between pa 10 3 and 
1.5 x 10 3 Z/q and its radius is pa 150 R Q , assuming in this case a distan ce of 


D = 1.8 kpc. The core mass is pa 0.50 — 0.55 M 0 (e.g., iBloeckerl 11993c Fig. 


[T]here). The companion will spiral-in due to tidal f orces when the primary 
radius is about quarter of the orbital separation (ISokerl. Il996l) . The initial 
orbital separation in this scenario is therefore ao ~ 3 AU. 

To avoid engulfment already on the RGB phase of the primary star, the 
primary size on the RGB must be < lO 2 i? 0 , which limits the initial mass 
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Figure 2: Light curve fit assuming a semidetached binary system made of a degenerate 
core and a main sequence companion, with the primary component filling the Roche lobe. 
B- and I-band light curves are shown in the lower panel using filled and open circles, 
respectively. Red lines are the best-fit models for each band. Residuals are shown in the 
upper panels. Note the different scale for each band. 


of the primary star to be 2.3 ^ M z amSi/Ms ~ 6.0 
Therefore, in this scenario the AGB star is more massive, and a more massive 
nebula is expected. 


fhben & Tutukovl . 119851 ). 


3.1.4■ Light curve 

Now the na tural question is whether this model is consistent with the 
observations of 1 Santander-Garcia, et al . (120151) . Clearly, in our model both 
components of the binary system have similar masses, and the stars are very 
close. Thus, the minima of the light curve will be broad, as a consequence of 
tidal distortion of the components of the bina ry system. This is in acc o rdanc e 
with the observed light curve obtained by ISantander-Garcfa et al. ( 12015! ). 
However, there is a fundamen tal difference between our model and that of 
Santander-Garcia et al. ( 2015 ). Due to the large luminosity contrast between 
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Primary 

Secondary 

M (M 0 ) 

0.495 

0.396 

R(Rq) 

0.506 

0.438 

r pole 

0.3748 

0.3226 

r point 

0.5230 

0.4044 

r side 

0.3949 

0.3406 

r back 

0.4250 

0.3669 


Table 1: Absolute masses and volumetric radii derived from the best fit to light curves 
assuming a semi-detached binary system. For each component, the radius relative to the 
semimajor axis is given in different directions accounting for deformation: towards the pole 
(fpoie) and towards the companion direction (r po i nt ), opposite (fback), and perpendicular 
to the orbital plane (r s id e )- 


the evolved object and the main sequence star, in our model the secondary 
star can be easily unnoticed in a first analysis, and the minima of the light 
curve are due to ellipsoidal variations of the post-AGB star. However, the 
light curve of this binary system can be easily confused with that of a binary 
system in which both components are post-AGB stars. 

To check whether our model is consistent with the observed data we 
computed synthetic light cu rves for this model. We used the latest ver sion of 
the Wilson-Devinney code ( Wilson fe Devinnevl . 1971 : Wilson! 1979 b This 
code allows to fit several light curves simultaneously providing parameters 
that are consistent with all the observed data. We have used the B- and 1- 
band light curves of Henize 2-428 to get the orbital and physical parameters 
of the system. Specifically, we obtained the inclination with respect the visual 
plane ( i ), the pseudo-potentials (Hi and fi 2 ), from which the radii of the 
components can be computed, the temperature of the secondary component 
(^efr, 2 ), the luminosity ratio for each bandpass (L 2 /Li), and a phase shift 
to account for period inaccuracie s. Sq ua re-root limb-darkening coefficients 
were interpolated in the tables of Claretl ( 2000a!) . and adjusted dynamically 
according to the current temperatures and surface gravities of the stars at 
each iteration. The reflection albedos were fixed at 1.0, appropriate for stars 
with radi ative en velopes, and the gravity darkening exponents were set to 1.0, 
following [Claret ( 200 1a). The best fit is obtained by differential corrections, 
iterating until the internal errors of the parameters are smaller than the 
corrections applied in three consecutive iterations and repeating this process 
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five times, selecting as best fit that with the smaller residuals. 

For this model we fixed the orbital separation to a = 1.27 Rq, which cor¬ 
responds to the observed period P = 0.1758 days. For the effective temper¬ 
ature of the primary we adopted = 45, 000 K, as explained in Sect. 13.11 
Finally, the mass ratio was kept fixed at q — 0.8. We fixed the primary com¬ 
ponent pseudo-potential to its critical value (Ox=3.42 according to the mass 
ratio) to force a semi-detached configuration and we tested fits with fixed 
values of H 2 hi steps of 0.1. This was done in order to prevent convergence 
towards an overcontact system if is set free. Nevertheless, the best fit 
was found close to the limit of overcontact with H 2 = 3.50. Figure [2] and 
Table Q] show, respectively, the light curves in the B- and I-band, and the 
relevant physical parameters of the best-fit solution. The best fit model has 
an inclination i = 69.792°, a phase shift of 0.0008, and the effective tem¬ 
perature of the secondary star is T e fj j2 = 41, 220 K. Finally, the ratio of the 
luminosity of the secondary star and the primary component is 0.683 in the 
B-band, and 0.679 in the I-band. Although the residuals of the fit are rather 
small (32.65 mmag in the B-band, and 23.51 mmag in the I-band), the top 
panel of Fig. [2] shows that there are systematic departures of the computed 
light curve from the observed one, which can be also easily noticed near the 
maxima and minima of the light curve in the bottom panel of this figure. 
We also tested different mass ratios, from 0.6 to 1.0, but the light curve fits 
were not significantly better. 

Clearly, the best-fit model is at odds with our initial hypothesis. Most 
significantly, the best-fit model shows that the luminosity of the secondary 
star is comparable to that of the primary component. This motivates us to 
look for other possibilities. 


3.2. An overcontact binary system 

Given that to fit the observational data the secondary component in 
the case of semidetached system should have a luminosity nearly equal to 
that of the primary we propose an alternative model in which both stars 
are at the same evoluti o nary stage. However, we relax the assumptions of 


Santander-Garcia et al .1 (120151) to allow for a more physically motivated so¬ 
lution. In particular, we assume that the position of Henize 2-428 in the 
Hertzsprung-Russell diagram corresponds to that of a model post-AGB star 
of mass ~ 0.5 M 0 , red line in Fig. [1] Thus, we adopt again T e ff, 1 = 45, 000 K. 
We note that since this mass is very similar to that found when a semi¬ 
detached system is considered we keep the orbital separation at a = 1.27 R e . 
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Figure 3: Light curve fits assuming an overcontact binary system made of two degenerate 
post-AGB cores. B- and I-band light curves are shown in the lower panel using filled and 
open circles, respectively. Red lines are the best-fit models for each band. Residuals are 
shown in the upper panels. Note the different scale for each band. 



Primary 

Secondary 

M(M e ) 

0.48 (0.22 - 1.50) 

0.47 (0.22 - 1.48) 

R (i? 0 ) 

0.57 (0.44 - 0.84) 

0.57 (0.44 - 0.83) 

r pole 

0.400 (0.386 - 0.420) 

0.397 (0.384 -0.415) 

r point 

contact 

contact 

r side 

0.431 (0.415 - 0.453) 

0.426 (0.412 - 0.448) 

r back 

0.495 (0.486 - 0.515) 

0.490 (0.486 -0.510) 


Table 2: Absolute masses and volumetric radii derived from the best fit to light curves 
assuming an overcontact binary system. For each component, the radius relative to the 
semimajor axis is given in different directions accounting for deformation: towards the pole 
(fpoie) and towards the companion direction (r po i nt ), opposite (rback), and perpendicular 
to the orbital plane (r s id e )- Numbers in parentheses indicate the range of possible values 
corresponding to fits with 0.8 < q < 1.2 and 1.0 < a < 1.9 Rq. 
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Additionally we did not take into account the mass ratio and the semi-major 
axis they derive because they are mainly constrained by the radial velocity 
fits, which we have argued that may not be related to the binary nature 
of the system. Nevertheless, the distance to the object and its luminosity 
may provide constraints on the scale of the binary system. Following our 
previous estimates in the case of a semi-detached binary, the total lumi¬ 
nosity of the system is scaled by L — 845(.D/1.4kpc) 2 L 0 . If we consider 
the independent distance measurements listed earlier in Sect. 12.21 Heniz e 2- 
428 may b e at d istance between 1.7 kpc and 2.7 ± 0.5 kpc (IMaciell . 11984 ; 


Frew et all 201511 . Then, the total luminosity of the system would be be¬ 
tween 1,246 and 4,415 L^. W e did not consider the distance estimate of 
( Santander-Garcfa et all 2015 1. which is shorter, to keep our estimates in¬ 
dependent of the light curve data we are using. Complementary, the total 
luminosity can also be computed from the effective temperatures, and relative 
radii derived from the light curve fits, along with the assumed semi-major 
axis. Therefore, the distance limits may provide constraints on the orbital 
semi-major axis. 

Taking all this into account we have performed fits to the light curves 
assuming different values of mass ratio (from 0.5 to 1.5) and semi-major axis 
(from 1.0 to 2.0 Rq), fixing T e s_\ = 45,000 K. The resulting light curves for 
our best fit model are shown in Fig. |3l and the physical parameters of the 
binary system are listed in Table [21 We note that, as expected, neither the 
semi-major axis nor the mass ratio have a significant impact on the quality 
of the fits at the level of precision of the light curves. However, they have 
an obvious impact on the absolute properties of the system. On the one 
hand, the absolute radius of each component is proportional to a, which 
defines the scale of the system. On the other hand, the relative size of each 
component (r, derived from the pseudo-potential) is correlated with the mass 
ratio because the stars of the binary are in contact. From the relative radii, 
the effective temperature and the semi-major axis, we computed the total 
luminosity of the system for each fit (L * (ar) 2 T e 4 ff ). The results show that 
the previously computed luminosity limits are consistent with a semi-major 
axis ranging from 1.0 to 1.9 _R 0 , almost independently of the value of the 
mass ratio. Given the period of the system, this can be translated into a 
total mass of the system in the range ~ 0.45 — 3.0 M 0 . If we a assume a 
mean measured distance of 2.2 kpc, then a ~ 1.3 Rq and the total mass 
would be ~ l.OMg, in agreement with the position of the system on the 
Hertzsprung-Russell diagram. 
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Given that the relative size of the components correlates with the mass 
ratio of the system, we have run fits by fixing a = 1.3 Rq and setting the 
mass ratio as a free parameter. The best fit is found for a = 0 .98 ± 0.20, 
very close to the ratio adopted by Santander-Garcfa et ah ( 20151 ). although 
it is poorly constrained by the light curves. As can be seen in Fig. 0 for this 
model the agreement between the theoretical light curves and the observed 
ones is excellent, both in the B-band and in the 1-band, being the residuals 
very small, 16.93 mmag and 3.25 mmag, respectively. For the best fit model, 
the inclination of the system is i — 63.2° and the effective temperature of 
the secondary star is T e ff i2 = 41,100 K. These parameters yield the masses 
of both components of the binary system, 0.48 M 0 , and 0.47 M 0 . Thus, the 
total mass of the binary system is 0.95 M 0 , well below the Chandrasekhar 
limiting mass. Hence, even if the system is made of two post-AGB stars, the 
total combined mass of th e system will not be en o ugh t o produce a Type 
la supernova outburst, as Santander-Garcfa et al. ( 2015f) claimed. In any 
case, stars with these properties are most probably post-AGB stars, and the 
component of lower mass is likely to be a helium white dwarf. 

Nevertheless, we stress here that the absolute properties of the compo¬ 
nents of Henize 2-428 are only roughly constrained by the light curve fits if 
we take into account the range of mass ratios and orbital semi-major axes 
that produce comparable fits to the light curves. Table [2] shows the properties 
of the best fit described above, with the range of possible values associated 
to the mass ratio and semimajor axis uncertainties. Although the combined 
analysis of light curves and the position in the Hertzsprung-Russell diagram 
favor a scenario of two similar stars with masses ~ 0.5 M 0 , only precise radial 
velocities can better constrain the model. 


4. Summary and conclusions 


We addressed the recent claim made by ISantander-Garcfa et al.l (120151) 


that the central binary system of the planetary nebula Henize 2-428 is a SN la 
progenitor in the frame of the double-degenerate scenario. This claim has at¬ 
tracted some attention, and if true it would be the first super-Chandrasekhar 
mass binary white dwarf system with short period discovered ever. More 
interestingly, it would be definitely located in the center of a planetary neb- 
ul a, prov idin g stron g evidence for SN la occurring inside PNe, termed SNIPs 
( Tsebrenko fe Soker . 2Q15h . However, we found that, although certainly pos¬ 
sible, this claim is not yet univocally supported by the observations, and that 
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different explanations are possible, and need to be considered. 

Onr statement is based not only on purely theoretical considerations 
about the evolutionary properties of the central object of Henize 2-428, but 
also on the way the observed data set has been analyzed. In particular, in 
Sect. [2] we argued that the claim for an equal-masses binary system and the 
required mass, luminosity, and radius of the two stars does not comply with 
evolutionary tracks of post-AGB stars (see Fig. [[]). We also analyzed criti¬ 
cally the way in which observations were interpreted, and we questioned the 
two most important assumptions made to study the observed properties of 
the binary system at the center of Henize 2-428. These are that both compo¬ 
nents of the binary system have exactly the same mass, and that, moreover, 
they are at the same evolutionary stage. 

In particular, in Sect. l2.2l we argued that there are strong th eoret ica l arg u¬ 


ments that pose a problem to the scenario proposed bv ISantander-Garcia et ah 
(120151 ). One of these arguments is that that the initial mass difference 
in the main sequence must be very small. The second one is that the 
population of twin binaries is small. Furthermore, in Sect. 12.31 we criti- 
call y exa mined the joint an alysis of the light curves and the spectrum of 


Santander-Garcfa et ahl (120151). S pecifically, we argued that the explanation 


of Santander-Garcia et al. (2 0151) of the variability of the He II 5412 A ab¬ 
sorption spectral line as arising from Doppler shifts of two absorption lines, 
one from each star, is not the only possible one. Instead, we suggested 
that the variability of the He II 5412 A spectral line can be accoimted for 
by a, possibly time-varying, broad absorption line from the central star on 
top of which there is a time-varying narrow emission lin e from the compact 
nebula or even from much closer to the star. Dobrincic et al.l ( 2008f) find 
the age of the Henize 2-428 equatorial ring to be 4,300(^0/1.8 kpc) yr. Al¬ 
though the ring is old, there is a compact dense nebula near the central star 
( Rodriguez et all 2001 ). It is quite possible that the He II 5412 A narrow 
emission spectral line sitting on top of the broad absorption line originates 
in the compact nebula, or from an outflow that feeds the compact nebula 
much closer to the star. This explanation could also be compatible with the 
observed spectrum, but more detailed studies and better observations are 
needed to resolve in detail the spectral feature. 

All these considerations led us to judge that although the first of the as¬ 


sumptions of Santander-Garcfa et ah ( 20151) - namely, that both stars have 
approximately the same mass - is quite possible, the second one - that both 
components of the binary system are at exactly the same evolutionary stage 
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- is not fully justified, and needs an independent evaluation. Since these two 
assumptions are quite restrictive, we explored alternative models. Accord¬ 
ingly, we analyzed other possible explanations, less extraordinary but that 
could fit equally well the observations and, simultaneously, do not contradict 
stellar evolutionary results (Sect. |3|. 

The first of these scenarios still hold the existence of a binary system, as 
the two (almost) symmetric broad minima in the light curve are attributed 
to tidal distortion caused by a companion of similar mass. In particular, we 
mentioned a binary system composed of either a post-RGB or a post-AGB 
star with a low-mass companion. In the first case, a low-mass main sequence 
star truncates the evolution of a star of initial mass M Z amSi ~ 1.0 M 0 on 
the upper RGB. In the second scenario, a main sequence star truncates the 
evolution of a star of initial mass MzamSi ~ 2.5 — 3.OM 0 on the lower AGB. 
However, we found that although this model is plausible it does not fit well 
the observed light curves. Consequently, we studie d a se cond possibility in 
which the assumptions of Santander-Garcfa et al.l ( 2015 1 are relaxed. We 
found that a binary system made of two post-AGB stars of masses 0.453 M 0 
and 0.437 M 0 can fit equally well the observed light curve of the system. 
Hence, if our interpretation of the observed data is correct the combined 
mass of the binary system would be below Chandrasekhar’s limiting mass, 
and it would not explode as a SN la. We conclude this short study by 
stating that the exact nature of the central object of Henize 2-428 is still to 
be determined. 
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